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A  New  Idea  for  Realizing  Positive  Rea! 
Immittance  Functions  of  Even  Rank 

KURT  H. HAASE 


Abstract 


A  novel  electrical  circuit  designed  according  to  the  theory  postulated  in  this 
paper  delivers  a  prescribed  immiltance  response-  A  driving-point  function  F(s) 
is  realized  by  a  single-T  circuit  and/or  ?.  single--  circuit,  implying  a  negative 
resistor  but  otherwise  positive  elements;  a  positive  real  immittance  function  F(s) 
is  realized  by  the  same  circuits  augmented  by  an  immittance  kf(s)  of  low  rank  at 
the  input. 

The  circuit  has  high  reliability  and  a  rugged  stability  in  the  face  of  environ¬ 
mental  changes.  In  addition,  its  low  weight  and  small  size  fulfill  requirements  for 
military  applications. 
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illustrations 


Terminated-T  Circuit  (a)  and  Term  mated -z  Circuit  (b)  Where 
F(c)  is  the  Driving-point  Immittanee 

Notation  (a)  for  Terminated-!  Circuit  (Impedance  Branches) 
and  (b)  Jor  Terminated-?  Circuit  (Admittance  Branches) 

ComparisoS  of  Networks  Having  the  Same  Driving-point 
Immittanee  in  Each  Column 

Two  Equivalent  Circuits  Realizing  the  Driving-point 
Impedance  F(s)  (v=4;  grounded  negative  resistance  R^J 

Example  2 

Two  Equivalent  Circuits  Realizing  the  Driving-Point 
Impedance  F(s),  With  (v=3).  of  Example  No.  3 


A  New  idea  for  Realizing  Positive  Real 
Itnirtiilance  Functions  of  Even  Rank 


i.  INTRODUCTION 


The  positive  real  (PR)  function  originally  defined  by  Oito  Brur.e  (1931) 
represents  a  general  class  of  rational  functions  composed  of  several  subclasses. 
The  subclass  treated  in  this  paper. 


F(s)  = 


N(s> 

D(s) 


,sv_1  +  ... +N.S+K 
_ v-l _ L _ 2 

sv  +  V  .sv~1  +  ...+D.s  +  D 
v- 1  1  o 


(1) 


is  a  function  of  even  rank  (2u),  and  denoted  as  a  function  of  the  ER  class.  All  its 
coefficients  N-  and  D-  are  nonnegative.  For  convenience  we  assume  that  the 
polynomials  N(s)  and  Bis)  are  normalized*  by  ■  +1. 

The  function 

N{s)  su+N  ,sv'l  +  ...  +  N1s  +  '5 

F(s)  r  —  =  J  -7ZI - .  (2) 

D(s)  sl+D  ,sl  1  +  . . .  +D.S  +  D 
v- 1  i  o 


(Received  lor  publication  4  December  1969) 

♦Unless  otherwise  stated,  the  term  polynomial  as  used  in  this  paper  always 
includes  normalization  in  this  sense. 
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considered  as  belonging  to  a  subclass  of  the  ER  class,  is  denoted  as  a  function  of 
the  H  class.  This  function  F  (s),  treated  in  Sec.  1,  has  some  properties  that  F(s) 
does  not  have. 

We  will  show  that  the  function  F(s>  can  be  realized  by  the  circuits  in  Figure  1,* 
where  it  is  the  driving-point  immittance.  The  terminating  immittances  Z(s)  and 
the  immittance  branches  X(s)  are  positive.  One  of  the  resistances  (conductances) 
Ru,  Rv,  Rw  is  negative;  It  has  to  be  realized  by  any  suitable  active  device.  If 
the  degree  v  is  an  even  one,  then  X(s)  and  Z(s)  are  PR  functions  of  rank  v- 1;  if 
the  degree  v  is  odd,  then  one  immittance  is  of  rank  i>  and  the  other  is  of  rank  v-2. 
We  shall  later  assume  that  F(s)  is  known  by  its  coefficients,  and  design  the  circuit. 

In  Sec.  2  we  will  show  that  it  is  possible  to  derive  an  K-class  function  F(s) 
from  an  ER-class  function  F(s)  of  the  same  rank  by  adding  an  impedance  of  low 
rank  in  series  or  in  shunt  to  the  circuit  that  realizes  F{s).  For  any  degree  v  the 
circuit  realizing  F(s)  must  be  a  single  T  or  single  terminated  at  both  ends. 

2.  REU.I/VHON  \M)  DEFINITION  OF  \N  H-CI.  \sS  FI  ACTION  Fi>) 

2.1  The  Norm  Function  F(s)  »ml  Its  Kcaiizifion 

Consider  the  circuits  in  Figure  2.  With  h  a  positive  constant  and  V  a  positive 


immittance,  let 

U  =  V(n-i)  (%) 

and 

W  =  -V(n-1  )/r  ,  (3b) 

Then 

1/U  +  1/V  +  1/W  =  UV  +  VW  +  UW  =  0  „  (4) 

XJ  +  V  =  nV  ,  (5a) 

V  +  W  =  V/n  .  f5b) 

U  +  W  =  V(o-l)2/n  .  (5c) 


*  In  this  figure,  as  in  all  others  in  this  report,  the  terminated-  T  circuit  (a) 
implies  impedance  branches,  ana  the  terminated-*  circuit  (b)  implies  admittance 
branches. 


fiw  driving-point  immittance  of  each  of  the  circuits  in  Figure  2  is 


.  ZOJ+V)  +  XfU+W)+XZ  _  VZfi2 + Y'Xin- 1 +  XZ5 
F*s>* - Jv+WT+XTZ -  - - V+xITzh 


«>) 


With  positive  constants  v,  x,  and  2.  normalized  PR  functions  n(s)  and  £>{s),  ant 
f  (s)  *  1  (dummy  function),  we  introduce 


V  =  vd(s) , 

ffa) 

X  3  xd>(s) , 

flo) 

Z  =  z|Cs)  =  z  , 

flc) 

set 

,  ._2 

z  »  1/n  . 

(7d) 

Then  by  substitution  in  Eq.  (6), 


F(s)  * 


dfe)-r3cCi-l)2dfe)^fe)r-xOfs)/v5 
9(sT+  l/vn+x54>(si/v 


(6) 


We  denote  the  form  of  the  function  presented  in  Eq.  (8)  as  the  norm.  Its  realiza¬ 
tion  according  to  Figure  2  needs  an  inverter  circuit  for  U  (which  is  negative  if 
n  <  1)  or  for  W  (which  is  negative  if  n  >  I).  In  both  events  the  inverter  presents 
the  technical  disadvantage  of  being  frequency-dependent  since  U  and  W  imply  cts). 

2.2  laiittMcr  Factions  F^(s)  ami  Fg{s*  to  be  Realized 

In  an  earlier  paper  (1966)  we  showed  that  F(s)  remains  invariant  when  in  any 
pair  of  the  functions  pis).  4>fs),  and  £(s),  the  functions  are  interchanged,  or  when 
the  functions  are  replaced  by  their  inverses  and  then  interchanged,  provided  that 
th;  original  constants  n.  v.  x  are  at  the  same  time  transformed  to  new  constants 
£',  v',  and  x!»  For  our  present  purpose  we  interchange  $(s)  and  f  (si  and  obtain 
F^fe)  *  F(s).  We  then  interchange  the  inverses  of  these  functions  and  obtain 
F^(s)  1  l/F(s).  The  transformation  formulas  are  listed  in  Table  1.  Wfcat  sc 
have  achieved  is  that  in  F^(s)  and  Fg(s),  V*  *  v';  concomitantly.  O'  and  W*  become 
mere  constants  and  thus  independent  of  frequency.  But  9(5)  and  <£>  (s)  become 
associated  with  the  constants  z’  and  x*,  which  are  cer vainly  positive.  Equivalent 
block  circuits  realizing  Ffe)  *  F^fs)  and  F(s)  *  Fg(s)  ar-i  presented  in  Figure  3. 
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Table  1.  Transformations  of  the  Constants  of  F(s)  for  Fajs) 
(<{>—*£,  <&)  and  l/FB(s)  (l/p  — 1/5,  l/<3» 


FA(s)  =  F(s) 

J/FB(s)  =  F(s) 

n1 

n/ (n-1) 

n/(n-l) 

u1 

l/n 

n 

V* 

(n-l)/n 

nffi-l) 

w' 

(l-n)/n2 

1-n 

x' 

x(n-l)/n2 

•/x 

z* 

v(n-l)2/n 

l/vn 

2.3  The  I’ll  uml  Normalized  Functions.  <4 ( s~)  and  <l‘(s) 

Let  a(s)  and  b(s)  be  two  polynomials  of  degree  ra>  and  let  o-(s)  and  fi(s)  be  two 
polynomials  of  degree  r  .  Assume  that 

^(s)  =  sb(Si£  PR  (9) 

and 

*(s)  =1^)  is  PR  .  (10) 

Substituting  these  functions  in  Eq.  (P.)  yields 


F(s)  = 


9 

s“a(s  Ms)  +  sodsjbfsj/vn  +  xnb{s)/i(s)/v 


(11) 


The  polynomials  in  the  numerator  and  denominator  of  Eq.  (11)  are  of  degree 
2  +i/  +r  and  normalized.  Their  degrees  are  the  same  as  t’->ose  of  N(s)  and  D(s) 

c\  li 

in  Eq.  (2)  when  we  provide  that 

for  even  degree  v,  v&  ~  vQ  -  —  (r  -  1)  ;  (12) 

for  odd  degi'ee  v,  (13a) 

a  6  *■ 

</„=!(</ -3). 


(13b) 


2.1  Derivation  of  the  Constants  n,  v,  \,  and  of  the  Functions  <$(s)  and 
<t>(s)  From  the  Coefficients  of  N(s)  and  H(s)  of  h(s) 

Assume  now  that  we  know  the  2v  coefficients  N  , . . . ,  N  ,,  and  D  , .... 

^  O  l/“i  w 

D  ,  .  We  want  to  find  the  three  constants  n,  v,  x,  and  the  v  +  v  coefficients 
v~  1  a  o 

of  £(s)  and  <Ms),  a  total  of  3  +  ^a  +  ra  unknowns.  Table  2  offers  4  +  v^  +  va  mean'-nS‘ 
ful  equations,  one  more  than  we  need.  This  surplus  equation  holds  only  when  the 
system  of  equations  is  consistent,  which  is  the  case  if  (a)  the  proposed  coeffi¬ 
cients  are  actually  those  of  an  H-class  function  F(s),  and  (b)  if  our  assumptions 
are  correct.  (We  shall  test  both  conditions  later,  j  Unfortunately,  as  a  glance 
at  v  -  4  shows,  the  equations  in  Table  2  are  nonlinear.  The  more  y  increases, 
the  more  complicated  they  become.  To  be  able  to  cope  with  any  degree  v,  we 
have  to  find  another  way  of  determining  the  unknowns. 

Note  that  for  any  degree  v, 


No  *  xboVvn  ' 

(17a) 

5o  =  xHboVv  • 

(17b) 

Therefore, 

ii  *  +  \Jdq/n o  . 

(18) 

Thus,  the  Constant  n  is  immediately  known. 

Let  us  define  the  polynomials: 

P(s)  *  sa(s)  -  b(s)/vn(n-l)  , 

(19) 

Q(s)  =  so(s)  +  xn(n-l)/3(s)  , 

(20) 

d(s)  *  (Nfs)  -  D(s)]/(n-l)  . 

(21) 

The  polynomials  P(s)  and  Q(s)  are  normalized,  but  cF(s)  is  not. 
and  Vq  are  defined  as  follows: 

The  degrees 

If  v  is  even,  t/p  °  =  £  v  • 

(22) 

If  v  is  odd,  Vp  =  j  (y+1);  Vq  =  • 

(23) 

The  degree  of  d(s)  is  v-1  , 

(24) 

$ 


& 

L: 

$ 


g  ! 


X  . 

fx  » 


j 

I 


f 

i 

i 


s 

i 

t 


tt 

i 

£ 

E 
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Since  the  coefficients  of  $(s)  and  <${s)  are  positive, 

P(s)  lias  only  positive  coefficients  if  n  <  1  ,  (25) 

Q(s)  has  only  positive  coefficients  if  n  >  1  .  (26) 

From  the  identity  of  Eqs.  (11)  and  (2)  it  follows  that 
P(s)Q(s)  =  |  nN  (r)  -  D(s)J/(H-l) 

»  sv  +  (PQ^.jS*'" 1  +  . . .  +  (PQ),s  4  PqQ0  .  (27) 

The  polynomial  product  P(s)Q(s)  is  known  in  its  summation  form.  By  solving  the 
equation  P(s)Q(s)  =  0,  we  can  transform  it  into  the  product  form.  In  this  form  it 
consists  of  the  product  of  some  linear  polynomial  factors  (corresponding  to  the 
real  roots)  and  some  quadratic  polynomial  factors,  each  of  the  latter  having  a 
negative  discriminant  (corresponding  to  the  conjugate  complex  roots).  From  this 
product  we  are  able  to  attribute  some  factors  to  P(s)  and  others  to  Q(s)  according 
to  the  information  given  by  Eqs.  (22),  (23),  (25),  and  (26).  But  some  ambiguity- 
may  remain  since  some  factors  may  be  either  in  P(s)  or  in  Qfs).  As  we  shall  s_e 
later,  this  ambiguity  can  be  cleared  by  using  the  surplus  equation. 

As  an  example,  assume  th»t  for  a  biquartic  function  F(s)  (i/=4),  with  n  >  1, 

P(s)Q(s)  =  s4  +  (PQ)3s3  +  (PQ),s2  +  (PQ),s  +  PoQo 

=  (s-s^)(s+sji(s+s,)(s+s3)  ,  (28) 

where  s^.  Sj,  s9,  and  s3  are  positive,  and  s£  <3,  <  Sg.  It  is  clear  that  (s-s^)  is 
a  factor  in  P(s)  only  when  n  >  I.  One  of  the  other  three  factors  is  also  in  P(s). 

We  thus  have  three  choices  for  distributing  these  factors  over  P(s)  and  Q(s),  as 
shown  in  Table  3. 


Table  3 


1  Choice 

Po 

P1 

^o 

No.  1 

"S4S1 

sr*4 

S2S3 

s2+s3 

No.  2 

~S4S2 

s2's4 

sls3 

sl+s3 

No.  3 

-®4S3 

S3-S4 

S1S2 

Sl+S2 

» 


m  i 


As  another  example  let  us  consider  a  bicubic  function  F(s)  (v=3>  for  which  we 
may  have  found  that 

Pis)Q(s)  =  s3  +  (PQ)0s2  *■  (PQ),s  +  P  Q 
-  I  o  o 

=  (s-s4>(s+Sj)(s+s2)  .  (29) 

where  s3*  Sj,  and  s2  are  positive,  and  Sj  <  s2-  In  this  example,  if  n  <  1,  then 
(s-Sj)  is  a  factor  in  Q{s)  and  the  other  two  factors  are  in  P(s).  But  if  n  >  1.  then 
we  are  left  with  the  two  choices  in  Table  4. 


Table  4 


Choice 

P  P, 

o  1 

% 

No.  1 

No.  2 

~sls3  s3'sl 

“S2S3  s2-s3 

s>  i 

Before  going  on  to  clear  these  ambiguities,  we  find  by  some  trivial  algebraic 
operations  that 

d(s)  =  (n+l)P(s)Q(s//n  -  s{a(s)P(s)/n  +  a(s)Q(s)J  , 
and  by  ordering  this  equation  we  get 

a{s)Q(s}  =  -o(s)P(s)/R  +  (n+l)P(s)Q(s>/{sn-d(s}J/s  .  (30) 

The  only  unknowns  in  Eq.  (30)  are  the  positive  coefficients  of  c(s)  and  o{s).  On 
both  sides  of  this  equation  there  are  normalized  polynomials  of  degree  r-1. 
Comparing  coefficients  yields  r- 1  meaningful  equations  by  which  we  are  able  to 
determine  the  v-2  unknowns.  Here  we  also  have  a  surplus  equation.  In  contrast 
to  the  results  of  the  earlier  coefficient  comparison,  however,  the  equations  de¬ 
rived  from  Eq.  (30)  are  linear  and  can  be  solved  bv  applying  Cramer’s  rule.  Or. 
the  right  side  of  Eq.  (30)  we  have  a  comparison  in  s' 1  that  has  no  match  on  the 
left  side.  It  yields  the  triviality 

(n+l)P0Q0  -  ndo  =  0  .  (31) 
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g* 
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When  all  the  unknown  coefficients  of  a(s)  and  o{s)  have  been  determined,  we  find 
the  remaining  unknowns  v  and  x  by 


l/vn(n-l)  =  a  ,  -  P 
v- 1  v 


xn(n-l)  =  Q  ~  a  , 
v  v  -i 
a  a 


{ 32 ) 


(33) 


The  coefficients  of  b(s)  and  $(s)  can  be  obtained  by  comparing  coefficients  in 
Eqs.  (19)  and  (20). 

2.3  Infractions  for  Sofvins  lie  Systea  cf  Linear  F^ualims 

We  assume  familiarity  with  Cramer’s  rule  (see  Hildebrand,  1956,  among 
Others).  The  derivation  of  the  matrix  system  to  which  the  rule  is  applied  deserves 
some  discussion.  These  matrixes  can  be  written  almost  immediately  since  their 
elements  have  to  be  taken  from  the  coefficients  of  P(s)  and  Q(s)  for  each  of  the 
possible  ambiguous  choices. 

Each  equation  derived  from  (30)  is  a  comparison  of  coefficients  associated 
with  sl  and  has  the  form 


j?o  +  k?o  (n+3HPQ)i+i/F'3i+i 


(34) 


where,  by  definition. 


Vj  =  0  for  ‘a  ’ 


a.  ,  =  0  for  i-k>  v  , 

l-K  a 


P^  *  0  for  k  <  0  , 

Q.  =0  for  i  <  0  . 
J 


(34a) 
(34b) 
(34c> 
(34  d) 


Suppose  we  write  these  equations  in  sequence,  starting  with  i  *  v-2  at  the  top 
and  ending  with  i  *  1  at  the  bottom  (we  need  only  these  equations  io  determine  the 
v  -2  unknowns).  Because  of  the  normalizations  there  are  some  constant  terms 
that  we  must  transpose  to  the  right  side,  obtaining 


* 


t 


12 


V1 


3 

where 


L  VjQi  +  .£  VkV“  *  Cj  .  A..  J  -  di+1 


nA.+  1  *  (t»+l>(PQ)i+I  -  &  +  nQi+i-t  >  •- 

a  a 

3Hr^i-5i+!>/(5-!)- 


(35) 


(35a) 


(35b) 


From  the  left  sid-:  of  these  equations  we  derive  a  square  matrix  §Q.  .,  P.  .j) 

i#  J  ^ 

of  v-2  rows  and  columns.  This  matrix  has  a  group  of  columns  listing  the  co- 
efficients  of  Q{ s)  and  another  group  listing  the  coefficients  of  P(s).  The  first 
column  in  each  group  lists  the  coefficients  sequentially  downward,  starting  at  the 
top  with  Qi-q  *  1  and  Pt,p  =>  1  .  The  list  is  repeated  from  one  column  to  the  next 
but  shifted  downward  by  one  position.  The  places  thus  vacated  then  contain  zero 
elements  - 

From  the  right  side  of  the  equations  we  derive  a  column  matrix  I»C  .§  of  i>-2 
rows.  The  difference  between  two  column  matrixes,  this  matrix  is  expressed  by 

ICj  *  ijA.+  Ig  ’  l3i+li  (36) 


With  these  column  matrixes  known,  we  can  apply  Cramer’s  rule  to  get: 


6  -det  |Qi  ..  P,  k|  ; 

*a.  'det  fci.y  Pi.k‘  * 


with  column  j  replaced  by  (3G); 


d«  pi.kl 


with  column  k  also  replaced  by  (36)  . 
Then 


(37a) 

(37b) 


(37c) 


ai  *  /A  ;  (38a) 

J  3 

and 

ffw/n  '  A  }&  -> 
s  ck 


(38b) 


2SLS 


r 


|i 
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2.6  The  Equation  of  Consist  C8e>  and  the  Definition  of  ll-Ciass  Functions  F(s) 

As  used  here,  the  surplus  equation  in  our  system  is  the  one  that  compared 
the  coefficients  associated  with  s°.  When  w_  substitute  the  coefficients  obtained 
through  Eqs.  (38)  we  obtain 


i>- 1  v~\ 

g,  Vi-  g,  3icix0* 


(39) 


for  which  A.  and  are  given  by  Eqs.  (35/, and  Cj  are  the  determinants  multiplied 
with  cofactors  of  a  square  matrix  that  is  obtained  When  the  square  matrix  of  v-2 
rows  and  columns  is  extended  by  the  row  i*0  and  one  row  is  deleted.  Assume,  for 
example,  the  degree  v-6.  The  extended  matrix  in  this  case  is  then 


1  i  iV 

i 

_£j 

1 

0 

k 

i 

i  t  , 

V  1 

Q3=1 

0 

P^l 

41 

0 

| 

£■  * 

£  ' 

3 

Q, 

i 

P2 

1 

(40)  | 

2 

*  ’ 

2 

^2 

P1 

P3 

\ 

1* 

J 

% 

Po 

P2 

r 

%  ! 

I  i 

0 

0 

Qo 

0 

P1 

1 

a 

! 

\ 

The  coefficients +Cj,  -c.,,  +  Cg.  - c^,  and  +  c-  are  obtained  by  respectively  deleting 
the  consecutive  rows  0,  1,  2,  3,  and  4.  The  coefficient  cq  is  defined  as 


‘'-.l 


co  = 


-A  -C./-1  , 

l  l  o 


(41) 


with  A  *  (r.*l)P  Q  'n  according  to  (35a).  Bv  (31)  we  can  add  (A  -  3_)c  on  the 

O  »)  O  --  Cl  o  o 


O  V)  O' 

left  side  in  (4  i)  and  obtain 


v-l 


6  =  -  52  cLc.  =  0  as  Ule  Equation  of  Consistency. 
i=0  1  1 


(42) 


Equation  (42)  yields  6=0  only  if  (a)  the  proposed  coefficients  N.  and  D.  are  those 
of  an  H-class  function  F(s),  and  (b)  the  correct  choice  has  been  made  in  selecting 
the  factors  of  P(s)  and  Q(s)  from  the  product  P(-.)Q( s)  in  Eq.  (27).  Thus,  this 
equation  also  defines  ihe  H-class  function. 
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Compact  formulas  for  the  coefficients  of  c.  can  b.'  given  for  the  relatively 
low  degrees  v  -  3  and  u  =  4. 

For  y  =  3, 

c2  =  -Qo'  ci=1*  co  =  -  f"Po  +  Qo(Pl  '  Qo)l/<*  +  l)Po%  '  <43) 

For  v  =  4  , 

c3  =  PlQo  '  W*  c2  =  Po'Qo'  “1  =Q1  *  ?1  * 

co  =  -  f nP0(c2  +  CjQo)  +  Qq(c2  +  ClP,)J /(n  +  1)P0Q0  .  (44) 

2.7  Nuaerical  Bxaaplc  No.  1  (Biquartic  H -Class  Function) 

Let  the  function  F(s)  of  degree  v  =4  have  the  coefficients  listed  in  Table  5, 
from  which  we  compute  the  coefficients  of  <T(s)  by  Eq.  (21).  With  n  *  2.  0.  we  use 
Eq.  (18)  to  obtain 

P(s)Q(s)  *  s4  +  15s3  +  63s2  +  41s  -  120 

=  (s-  l)(s+3)(s4-5)(s+8)  .  (45) 


Table  5 


i 

Ni 

D. 

i 

3. 

1 

0 

60 

240 

-180 

1 

118 

195 

-77 

2 

75.  5 

88 

-12.5 

3 

15. 5 

IS 

5 

4 

1.0 

1.0 

0 

The  three  choices  for  selecting  P(s)  and  Q(s)  from  this  product  are  listed  in 
Table  6  together  with  the  coefficients  c.  computed  by  (44)  and  6  by  (42)  for  each 
of  the  choices.  According  to  the  last  column  in  Table  6,  the  correct  choice  is 
No.  2,  for  which  by  Eqs.  (37)  we  obtain 

,  A  =  -45.5,  £  *  -49  . 

tt.  U 

0  o 


£  -»  -  7 


(46) 


Table  6 


Choice 

No.  1 
No.  2 
No.  3 


*0 

Qi 

c3 

c2 

C1 

Co 

6 

40 

13 

119 

-43 

11 

-4 

-351 

24 

11 

151 

-29 

7 

-i.4 

0 

15 

8 

169 

-23 

1 

0 

-128 

Then,  by  Eqs-  (38),  we  find 


a  =6. 5.  a  =  7  ; 
o  o 


and  by  (32)  and  (33), 
v  =  0.  2  .  x  =  2  . 

Coefficient  comparison  in  (19)  and  (20)  yields 


bo*2,  0o  =  e. 


The  two  equivalent  circuits  realizing  the  impedance  F(s)  are-  shown  in  Figure  4. 
Their  elements  computed  according  to  Tables  1  and  Figure  3  are  listed  in  Table  7, 


Table  7.  Circuit  Elements  in  Figure  4  for  Examples  Nos.  1  and  2 


Example  No.  1 


Example  No.  2 


FA(s)  =  F(s) 

l/FB(s)  =  F(s)  FA(sl  =  F(s) 

l/FgCs)  =  F(s) 

9. 5 

0.5  0.2745587 

0. 2745587 

0-  5 

0.  5  0. 72544 12 

0. 103S12S 

-0.  25 

-1.0  -0.1991762 

-0. 3784714 

2.0 

2, 0  2.  3095979 

0. 3308403 

169/180 

338/90  1. 6368449 

3. 1 103057 

0. 325  ; 

1.  3  0.  3733889 

0. 7095074 

13/90 

52/90  0.2708061 

0.5145811 

0.5 

j  0.5  0.3654515 

2.5513119 

0.  5/6 

|  0-5/6  0.0799116 

0. 5578634 

KfllM'1*  l""?!  *  U!  IWflf'tf  Jfc&JiAft'  lltc?  CttK/  < 


18 


I 


It  can  easily  be  shown  that  cT{ s),  defined  in  (2  0,  is 

cf(s)  =  (N(s)  -  D(s)|/(n -  0(1  - k)  .  (53) 

But  then  Eq.  (42)  can  be  expressed  as 


5  = 


v-1 

T,  c,d.  =  0  , 
i=0  1 


(54) 


where  this  time  the  coefficients  Cj  are  derived  from  (52)  instead  of  from  (27),  and 
n  is  derived  from  (50).  Since 

l/k  =  n“  [n"(l  -2/nQ)+  1 J / (n^  -n2)  ,  (55) 

the  constant  k  can  accordingly  be  replaced  by  n^.  Instead  of  solving  directly  for 
k  it  is  more  convenient  to  first  determine  nQ.  The  H-class  function  F(s)  in 
Eq.  (48)  has  Leen  found  when,  after  a  search  with  trial  parameters  nQ,  Eq.  (54) 
yields  6  =  0  We  show  this  in  Sec.  3.2. 

3.2  Numerical  Kxumplc  No.  2(ltiquarlic  H'CIass  Function) 

Let  a  biquartic  ER-class  impedance  function  be  described  by  the  coefficients 
listed  in  Table  8.  The  table  also  lists  the  coefficients  of 

d(s)  =  (N(s)  -  D(s)]/{»-  1).  (56) 

2 

The  function  F(s)  is  PR.  The  minimum  of  ReF(jt;  )  is  F(0)  =  l/n  =  0.4.  We  dis¬ 
regard  that  F(s)  -  F(0)  can  be  realized  by  the  well-known  ladder  procedure.  But 
we  want  to  point  out  that  any  positive  constant  0  <  k  <  0.  4  can  be  subtracted  from 
F(s),  and  so  the  value  nQ  can  range  between  2.  5  and  \ZTi  =  1.  5811388. 

Let  us  first  test  whether  F(s)  is  an  H-class  function.  (This  can  be  done 
either  by  subtracting  the  constant  k  =  0  or  letting  nQ  =  V2-  )  If  we  treat  F(s) 

as  though  it  is  an  H-class  function,  Eq.  (55)  will  yield  the  solutions 

s.  =  1.  2143093,  s,  *  3.2212775,  s„  =  4.7705955,  s,  =  8.  1341325.  (57) 

With  these  solutions  we  obtain  the  results  listed  in  Table  9.  As  the  last  column 
in  Table  9  shows,  F(s)  is  not  an  H-class  function.  We  therefore  use  a  series  of 
trial  parameters  nQ  to  compute  the  6j,  6^,  6 ^  for  the  possible  choices.  Some  of 
the  results  are  presented  ir.  Table  10.  As  the  table  shows,  choice  No.  2  is  the 


f 

] 


i 


Table  8 


-247.7893405 
-105.9987734 
-17.2075930 
-0. 6883037 
0 


Table  9 


Choice 

Po 

P1 

Qo 

«1 

No.  1 

No.  2 

No.  3 

-3.9116272 

-5.7929785 

-9. 8773527 

2.0069682 

3.  5562862 

6.9198232 

18. 8046559 

26.2022980 

15.3674119 

12. 9047280 

1L  3554100 

7.9918730 

Table  10 


1.6 

1 

1.  66 

1.68 

1.88 

1.  89 

-409.7952392 

-476.4259531 

-496. 3981346 

-658. 3761346 

-665. 1590852 

-8.4908898 

-0. 5233455 

+0. 8875778 

+0. 3571147 

-0.  0203837 

-195.0626017 

-175.4795503 

-169.4790472 

-17.0.  4249847 

-138.3943812 

f 

correct  one.  The  value  6 2  changes  its  sign  of  polarity  from  nQ  =  1. 66  to  1.  68 
and  again  from  1.  88  to  1.  89.  The  exact  zero  crossings  of  62  can  be  found  by  any 
well-known  interpolation  formula  (we  suggest  the  one  by  Aitken  in  Abramovitz  and 
Stegun,  1964).  The  results  of  interpolation,  together  with  the  constants  Y,  n,  and 
n,  are  presented  in  Table  11. 

For  K  =  0.  2  and  k  =  0.  8  the  function  F(s)  is  the  same  as  in  Example  No.  1. 

The  impedance  function  has  the  same  equivalent  realizing  circuits  shown  in 

Figure  4.  But  since  we  want  to  realize  KF(s),  the  resistances,  inductances,  and 

inverse  capacitances  listed  in  Table  7  have  to  be  multiplied  by  K  *  0. 2.  Finally, 

to  obtain  the  impedance  F(s),  both  circuits  have  to  be  augmented  with  a  series 

resistance  R  =  k  =  1-K  -  0.  8. 
o 


Table  11 


no 

15/9 

1.  8894704 

62 

-0. 0000013 

+0. 0000020 

n 

1.  5611388 

1.5811388 

n 

2.0 

3. 6422085 

K 

0.  2 

0.  3510830 

_ i 

With  K  =  0.  3510830,  the  constants  obtained  for  F(s)  are: 

n  =  3.  6422085,  v  =  0.0818905,  x  =  0.3919552  ,  (58a) 

and  the  coefficients 

aQ  =  6.0443381,  aQ  =  (i.  0029088,  bQ  =  2.  5409135,  pQ  =  5.  18 1825.  (58b) 

The  impedance  realizations  of  F(s)  are  again  those  shown  in  Figure  4,  with 
the  circuit  elements  listed  in  Table  7.  The  element  values  have  accordingly  to 
be  multiplied  by  the  factor  K  and  the  circuits  have  to  be  augmented  with  a  series 
resistance  of  Rq  =  1-K  =  0.  6489170  at  the  input  to  obtain  realization  of  the  im¬ 
pedance  F(s).  We  have  thus  found  four  equivalent  circuits. 

We  want  to  point  out  that  if  f(s)  •<  1,  the  constant  k  does  not  necessarily  have 
to  be  positive.  Since  we  have  accepted  one  negative  resistance  in  the  realization, 
we  can  as  well  accept  a  second.  There  is  also  the  chance  that  this  second  nega¬ 
tive  resistance  may  be  cancelled  out  by  a  positive  one  if  F(s)  is  part  of  a  larger 
circuit. 

3.3  The  Function  F(s)  Obtained  by  Subtracting  kf(s)  From  F(s) 

if  trials  with  an  assumed  function  f(s)  *  1  fail  to  yield  an  H-class  function,  we 
have  to  find  a  more  costly  function  f(s)  =  n(s)/d(s).  Considering  th?  decompositic  n 
(Eq.  (44)]  only,  let  d(s)  be  a  part  of  D(s), 

D(s)  =  d(s)D'(s)  .  (5.9) 

Assume  that  the  degree  of  n(s)  is  one  order  less  than  the  degree  of  d(s).  Then 
according  to  (44), 


KF(s)  =  (1-  k)N(s)/D(s)  , 


(60) 


IF  » 


r 


Ef 


2  i 


with 


N(s)  =  (N(s)  -  kn(s)D'(s)]/(l-k);  (60a) 

D(s)  ■=  D(s)  =  d(s)D'(s)  .  (60b) 

We  can  alsc  allow  n(s)  i o  be  of  the  same  degree  as  d(s),  changing  (60)  and  (60a) 
accordingly.  We  always  have  to  ensure  that  KF(s)  i.i  (44)  is  PR,  which  is  true 
when 


Re  F(ju)  -  Re  kf(jw)  >  0  . 


(61) 


f. 

* 

i  . 

f-  ■ 


fe  ' 


Let  us  now  consider  the  following  example. 

3.  t  Numerical  Example  No.  3  (Ilieuliic  H-Class  Function) 

For  the  coefficients  of  a  bicubic  impedance  function  F(s)  listed  in  Table  12, 
we  find  that  the  denominator  of  F(s)  is  j 

1 

D(s)  =  (s  +  so)(s2  +  D{s  +  D^)  ,  (62)  I 

where  ! 

=  0.  5207709,  D\  =  0*4792291,  D*  =  3.6004314  .  (62a)  I 

o  i  o  -  } 


I* 


i- 


Table  12 


We  have  chosen  f(s>  =  l/(s  +  s  ).  In  Figure  5  we  compare  Re  F(ju)  and  Re  f(ju) 

2  0 

versus  fi=  u  .  Since  F(0)  =  0.  928,  and  f{0)  =  1.  9202301,  the  constant  k  must  be 
less  than  9. 483.  Factorization  of  P(s)Q(s)  allows  two  choices,  6  j  and  6^,  whose 
values  for  a  series  of  parameters  k  are  listed  in  Table  13. 


23 


Table  13 


k 

0.48 

0.40 

0.  399 

0.  3995 

*1 

-0.2067047 

-0. 0007434 

+0.0115423 

+0. 0053377 

*2 

-0.5809142 

-4.0111988 

-4. 04S2474 

i 

-4-0287219 

1 _ - _ 1 

Being  satisfied  with  the  small  6  j  =  -0.  0007434  for  k  =  0. 40  we  go  on  to  obtain 
No  =  0. 3.  N1  =  8. 0,  N2  =  7.55,  N,  *  1. 0  .  (63) 

The  constants  for  F(s)  are 

n  =  2.  5,  v  =  2.  0,  x  =  3. 0  ,  (64) 

The  coefficients  of  $(s)  *  s(s+ao)/(s  +  bo>  are 

aQ  =  0. 8,  bQ  *  0. 5;  4>{s)  *  l/s  .  (65) 

The  two  equivalent  circuits  realizing  the  driving-point  impedance  F(s)  are  shown 
in  Figure  6,  and  the  element  values  are  listed  in  Table  14. 


Table  14,  Example  No.  3 


T  Circuit 

*  Circuit 

Ro 

0. 7680920 

0. 7680920 

R1 

0.24 

0. 24 

R2 

0.  36 

0.  16 

R3 

-0.  144 

-0.4 

R4 

2.304 

6.4 

L 

1 

1.725 

4.8 

4 

2.88 

8.0 

Co 

2.5 

2.5 

C1 

1/4. 05 

1/1.8 
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t.  COWIKNTS 

It  must  by  now  be  apparent  that  the  driving-point  realizations  discussed  in  this 
paper  have  some  similarity  to  the  well-known  Brune  realizations  (1931).  Our 
realization,  however,  is  a  single  terminated  two-port,  whereas  Brune’s  is  a  cascade 
of  several  two-ports.  The  single  two-port  allows  us  to  induce  the  negative  resist¬ 
ance.  Where  application  of  the  Brune  cycle  requires  that  a  duplex  zero  appear 
for  the  real  component  Re  F(jo),  our  method  requires  that  either  Eq.  (42)  or  (54) 
yield  6  =  0.  The  Brune  realization  allows  transforming  the  circuit  into  a  Bott- 
Duffin  (1949)  circuit.  We  have  not  yet  found  a  similar  equivalent  for  our  circuits 
even  at  the  expense  of  using  more  elements  to  avoid  the  negative  resistance.  We 
can  get  a  Bott-Duffin  equivalent  only  when  d(s)  =  l/<l>(s>;  in  such  case,  however, 

F(s)  is  nothing  more  than  a  frequency-transformed  biquadratic  function.  Although 
this  event  is  almost  trivial,  cur  realization  procedure  at  least  offers  the  means 
of  discovering  the  transformation.  Whether  an  H-class  function  can  always  be 
derived  from  an  ER-class  function  by  splitting  off  a  proper  function  kf(s)  will  have 
to  be  determined  empirically  through  a  computer-aided  design. 
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